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Abstract 

Extending  the  theory  vhich  developed  by  Kohn,  Onffroy  and  Rehr  for  a one- 
diiaensional  non-periodic  systec’.,  we  have  developed  a theory  that  describes  hov?  the 
electronic  density  distribution  near  an  insulator- insulator  interface  approaches 
the  bulk  behavior.  The  exponential  decay  constant  depends  on  the  band  struc- 
tures of  the  two  insulators  and  the  positions  of  localized  surface  state  energies. 
In  addition,  we  also  study  the  behavior  of  the  electron  density  perturbation  due 
to  a defect  across  several  successive  interfaces. 


* Svqjported  in  part  by  the  Office  of  Naval  Research  and  the  National  Science 
Foundation. 
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I.  Intxoduction 


'l^!r.:  i-  ’i.,  4d 


/ 


In  an  infinite  and  perfectly  periodic  crystal,  the  eigenfunctions  of  the 
one-electron  Schrodinger  equation  are  the  Bloch  vgaves.  The  electron  density 
distribution  is  a periodic  function  with  the  same  periodicity  as  the  potential. 

In  a crystal  with  an  interface,  the  periodicity  of  the  potential  is  lost.  The 
eigenfunctions  of  this  crystal  are  no  longer  the  Bloch  waves . In  the  vicinity 
of  the  interface,  the  electron  density  distribution  is  perturbed  by  the  exis- 
tence of  the  interface.  However,  de^  inside  each  of  the  crystals,  the  electron 
density  distribution  approaches  that  of  the  infinite  perfect  crystal.  In  this 
part,  we  want  to  find  in  \diat  manner  the  electron  density  distribution  approaches 

that  of  the  infinite  perfect  crystal,  as  we  move  away  from  the  interface  into 

1 2 

either  of  the  two  crystals.  Kohn  and  Onffroy  Rehr,  and  Kohn  developed  a 
procedure  to  construct  generalized  Wannier  functions  for  a system  with  an 
inpjrity  or  a surface.  They  found  that  the  generalized  Wannier  functions 
approach  the  Wannier  functions  of  the  perfect  crystal  exponentially  with  dis- 
tance from  the  defect  or  the  surface.  We  want  to  extend  their  idea  to  study 
the  perturbation  of  the  electron  density  distribution  by  an  interface  in  one- 
dimensional  insulators.  For  an  Insulator  the  Fermi  level  lies  at  an  energy 
corresponding  to  a gap  for  each  of  the  two  crystals,  i.e.  all  the  energy  bai 
of  the  perfect  crystals  are  either  full  or  enpty.  We  find  that  the  perturbed 
electron  density  distribution  approaches  that  of  the  perfect  crystal  in  an 
exponential  manner.  The  decay  constant  in  each  crystal  is  determined  by  the 
relative  positions  of  the  following  quantities:  the  positions  of  the  band  edges 
of  the  last  filled  and  first  unfilled  energy  bands  of  the  other  crystal;  the 
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positions  of  the  localized  surface  state  energies;  and  the  position  of  the 
branch  point  energy  \diich  are  located  between  the  last  filled  and  first  vmfilled 
energy  bands  of  this  crystal. 

We  also  extend  this  theory  to  study  more  general  cases,  the  vay  in  which 
the  density  perturbation  due  to  a defect  appears  across  several  successive 
interfaces.  Again,  we  find  that  this  perturbation  decays  exponentially  and 
obtain  a cocplete  prescription  for  determining  the  decay  constants  for  an 
arbitrary  succession  of  periodic  lattices.  We  believe  that  our  results  have 
some  promise  for  developments  of  a theory  for  the  following  interesting  chemi- 
cal polymer  problem:  if  we  replace  one  atom  by  a different  atom  at  one  point, 
how  does  the  perturbation  behave  along  the  chain  of  the  polymer? 

In  Section  II  we  will  derive  the  general  theory  of  the  electron  density 
perturbation  due  to  the  presence  of  the  interface  in  a one-dimensional  insulator 
of  two  lattices  in  contact. 

In  Section  III  we  will  use  this  theory  to  discuss  the  electron  density 
perturbation  due  to  an  isolated  defect  in  more  complex  insulators  having  n 
arbitrary  different  lattices  in  contact. 

Section  IV  contains  a sunnary  of  our  results. 


II.  Isolated  Defect  in  a Perfect  Insulator  and  Simple  Interface 
A.  General  Theory 

Before  proceeding  to  derive  the  theory,  let  \is  briefly  review  some  proper- 
ties  of  the  perfect  periodic  lattice.  The  eigenfunctions  of  the  one-electron 
Schrddinger  equation  are  the  bloch  waves  i|)^(x,k)  with  energy  E(n,k)  (eigen- 
values). Here  k is  the  irave  number  and  n is  the  band  index.  If  we  treat 
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k as  a ccxiplex  variable^  k * g + ih,  than  '1'^  (x,k)  and  E (n,k)  represerit 
the  branches  of  nultlvalued  functions  of  k . They  are  analytic  everywhere 
except  at  a set  of  branch  points  away  from  the  real  axis.  If  g is  restric- 
ted to  the  fundamental  interval  - — <g  < — , these  branch  points  have  the 
ir  n b ■ b 

form  k”  = — + ih  , where  (x^  * 0 or  1 for  n odd  or  even,  respectively. 

They  can  be  determined  fran  the  Kramer's  plot  (Fig.  1)  of  cos(l;a)  versus  real  E 
The  branch  points  k^i  correspond  to  the  (n  + l)th  extrenun  of  this  plot.  The 
energy  E^  = ECW^")  at  these  branch  points  are  real. 

We  now  consider  an  electron  in  a one-dimensional  potential  vMch  is  taken 
to  be  extended  over  the  interval  - d/2  < x < d/2,  with  Hamiltonian  (Fig.  2) 


H = - ^ + V (x) 
dX^ 

/v,(x^ 


V(x)  (x) 
\v.,(x') 


-d/2  < X <-N 


-N  < x<0 


0 < x < d/2 


where  v^(x>  and  V2(x)  are  potentials  with  the  periodic  and  syninetry  properties 


v^(x+b)  = v^(x); 


v^(-x)  = v^(x),  where  v=  1,  2. 


Here  b is  the  period  of  the  potentials  v^(x)  and  V2(x) , and  vhere  N is 
the  length  of  the  section  at  vMch  the  potential  is  Vj^(x) . 


The  eigenfunctions  <))  (x)  are  required  to  satisfy  the  differential  equation: 


H (j)(x)  = E (^(x) 


and  the  conventional  periodic  boundary  conditions: 


(}.(d/2)  = K-d/2)  ; 


4. '(d/2)  = (t.'(-d/2). 


viiere  <)> ' (x)  = (x) 

dx 


When  N equals  zero  (i.e.  the  case  of  the  perfect  crystal)  these  eigenfunc- 
tions are  the  Bloch  waves 

ip  (x,k  ),  where  k2  = — (-  — , — .O—.Cjb  (2.5 

2 2 ^ d 2b 

The  corresponding  eigenvalues  are  groi;q)ed  in  quasi-continuous  energy  bands,  each 
containing  d/b  levels . 

When  N is  not  equal  to  zero,  the  energy  states  of  the  original  band  split 
into  two  groiqjs;  one  having  m localized  impurity  bands  and  the  other  having 
(d/b-m)  minibands  (see  Fig.  3).  When  d these  localized  impurity  bands 

approach  the  bound  states  and  the  minibands  become  scattering  states  vdiose 

energy  gaps  approach  zero^  (see  Fig.  4)  . nrccFss'n»  U 
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then  we  know  that,  in  this  form,  (x,k2)  is  an  analytic  function  of  complex  k2 
except  at  the  branch  points  vAiich  connect  the  miniband  m to  its  neighbors 


nri-  1.' 


M 


n (x.k  ) =2  (x,k  ) , 


(2.11) 


m=0 


where  m runs  over  all  minibands  including  the  inpurity  bands  and  n (x,k2)  is 
analytic  at  all  branch  points  connecting  the  minibands  m^. 

The  total  electron  density  from  all  the  states  of  this  band  is 


n (x)  = 


i M ^ 

j dk^n  (x,kO  '’‘■'"z’ 


(2.12) 


as  d -*■  a>,  we  have 


in=0  _ TT 

' d 


TT 

n (jO  = (x,lc2j) 


(2.13) 


TT 

b 


0^  * 

where  n the  electrcxi  density  of  the  scattering  states,  is 


n®""  (x,k  ) = \ (X)  <1^  (X)  , (2.14) 

2 ^2  ■‘^2 

vhere  (x)  are  the  scattering  states,  (+)  and  (-)  denoting  states  with 

outgoing  and  incoming  scattered  states,  respectively. 

r?  (x,k2^),the  electron  density  of  the  bound  states,  is 
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(2.15) 


(x,k2^>=  (x,k2^)  (x,k2^^ 

where  (j)  (x,k2‘^)are  the  normalized  bourvd  state  wave  functions  with  wave  vector 
k2-^  whose  arrplitudes  are  determined  by  tlie  following  asymptotic  condition: 


(x,k2-^)=  \l)^  (x.kj^)  , x>o, 


(2.16) 


where 


is  determined  by  the  relation: 


B 


(2.17) 


and  where  eT  is  the  bound  state  energy  which  is  in  the  energy  gap  of  a perfect 
lattice  with  potential  (x) . Here  extends  over  all  the  occupied  bound  states . 

The  scattering  process  due  to  the  existence  of  the  potential  Vj^  (x)  can 
be  describee'  transfer  matrix  which  relates  the  amplitude  and  phase 

of  the  sc  ■'ves  to  those  of  the  incident  waves.  Here  the  superscript 

(2)  den  ^ as  an  independent  variable.  We  use  two  different  vjays  to 

derive  the  transfer  matrix:  (1)  Saxon's  nulti-scattering  method;'^  (2)  Heine's 

5 

Bloch  wave  matching  method. 

(1)  Saxon's  multi- scattering  method: 

We  can  regard  the  presence  of  the  potential  v^^  (x)  as  a succession  of 
scattering  centers  vdiich  are  located  from  x = -M  to  x = 0 (see  Fig.  2) . When 
N = 1 , the  solution  can  be  viewed  as  an  unperturbed  Bloch  vrave  interacting 
with  an  isolated  localized  inpjrity.  This  scattering  process  can  be  described 
by  a scattering  matrix  vMch  relates  the  amplitude  and  phase  of  the 

scattered  waves  to  those  of  the  incident  waves.  For  a finite  N , the  passage 
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of  an  electxon  throujih  the  N scattering  centers  can  be  thought  of  as  a nulti- 
ple  scattering  of  Bloch  electron  waves  by  those  scattering  centers.  The  trans- 
fer natrix  is  then 


M' 


(2) 


- R®)  s.i<»  (S.2<»)-1  R«)  - - - R<2>  S.(^_„«>(S.,/2>)-1  R<2>  S., 


(2.18) 


where  S is  the  translation  matrix  for  translation  to  the  right  through  a 
a 

lattice  constant:  (From  now  on,  we  assume  that  the  lattice  constant  b equals  one.) 


\diere 


has  the  following  properties: 


(2.19) 


^1 


(2) 


(V  = ^2/^^ 

Rii<2>  (V  = ^22^^^ 


(2.20) 


and  vdiere 


g(2)  , s^(2)  ^(2) 


(2.21) 


From  Appendix  A,  we  have 


(jO  = <p2  'i>2 

'('ll(x)  = a^  K^3^(x,k^)  + b^  ^2 

’^III  ~ ^3  '^2  ^3  ^2  (2.24) 

Using  tlie  continuity  conditions  at  the  interfaces,  we  obtain 

m(2)  L 

W ~^)  ) W (4)2(0>k2)  >'1^2^^’ ^ 

\ -521®’  5U®  / , e-«^-¥%2(»j(2.25) 

and 

(k2)  = (k^)  (2.26) 

vtere  w(f  ,g)  is  the  Wronskian  of  the  functions  f and  g : 

w (f,g)  =fg'  - f’g  (2-27) 

and  idiere 
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(2.28) 


= w (i|;2(0.k2).  ) 


=w  (i/;^(0,k^),i|)2(0.k2)  ) 


Q22^^^  =w  (ip^(0.k^),.^2<°*-^>  > 


For  a finite  N , the  ratio  of  the  transfer  matrix  elements  ^ 

analytic  function  of  1^2,  except  at  branch  points  v^ch  are  the  same  as  those 

(2)  (O') 

of  E'  (k2)  and  for  the  possible  pole  sirgularities  at  the  zeroes  of  1^2'  ^ • 

/0\  / '\ 

As  N goes  to  infinity,  though  the  retains  the  just  mentioned  branch 

point  structure,  two  additional  analytic  features  arise:  (1)  two  new  sets  of 
branch  points  occur  at  valvies  of  ^2  corresponding  to  the  band  edges  of  a per- 
fect lattice  with  potential  (x)  and  (2)  an  infinite  sequence  of  poles  arise 

(2) 

from  the  infinitely  many  zero-r-s  of  F^2'  ' • 

The  perturbation  of  electron  density  at  X2  (x2>0)  due  to  the  existence 
of  the  potential  (x)  is  (Appendix  B) 

/•H  M (2) 

6n(x2)  = / dk2  — - — (x2,k2'^)  (2.29) 

k ^2® 


We  would  like  to  discuss  Eq.  (2.29)  in  the  following  two  sinple  cases:  an 
isolated  defect  in  a one-dimensional  perfect  insulator,  and  an  interface. 
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B.  Isolated  Defect 


When  N = 1 , the  systan  can  be  thought  of  as  an  isolated  defect  at  the 
origin  in  a perfect  lattice  vd.th  the  potential  ^2  (^) . The  strength  of  the 
defect  potential  is  (x)  (see  Fig.  5). 

Since  regular  in  the  upper  half  of  tlie  plane  in  the  region 

betvjeen  the  real  axis  and  the  branch  points  . except  for  the  possible  pole 

(2) 

singularities  at  zeroes  of  1^2'  ' corresponding  to  the  existence  of  the  bound 

states,  we  can  free  to  shift  the  path  of  integration  of  Eq.  (2.29),  from  the 

real  axis  of  to  the  contour  ^ as  indicated  in  Fig.  6.  This  results 

because  of  the  periodicity  of  the  integrand  in  Eq.  (2.29)  (the  vertical  portions 
2 

cancel  each  other) . Thus 


^2^^  j 


-2Tri  ^ ^ res 

• I 

J 

where  j runs  over  all  the  occupied  bound  states  and  where  j runs  over  all 

the  poles  which  are  located  between  the  real  axis  and  the  contour  C . 

From  Eq.  (2.30),  6n  (X2)  will  depend  upon  the  strength  of  v^  (x)  and 

(2) 

V2  (x)  , the  positions  of  tlie  zeroes  of  M,2^  ^ and  the  positions  of  the 
(2) 

braxich  lints,  E^'  , which  are  located  between  the  last  filled  and  the  first 

unfilled  energy  bands  of  a lattice  with  the  potential  ^2  (jO-  We  can  classify 
all  situations  into  three  cases: 


Mj_2 


(2) 


1^0 


(2) 


(X2,k2^  )\p2  (^2’  ^ 


(2.30) 


Case  1:  E lies  between  the  energies  of  the  first  ijnoccupied  state 
o « 

(2) 

and  the  last  occi^ied  state:  the  residues  due  to  the  zeroes  of  ^^2 

2 

cisely  cancel  the  bound  state  contributions.  We  obtain 


C 

“ / 7^  '^2  "^2  V 


(2.31) 


This  leads  to  the  asymptotic  behavior: 


|6n  (^2)  1< (const.)  e~^^  ^ as  Xj 


(2.32) 


where  h ° is  the  distance  from  the  real  axis  to  the  branch  points  k ° in 
2 2 

the  k2  plane. 

Case  2:  lies  below  the  energy  of  the  last  occupied  state:  the 

residues  from  the  zeroes  of  ^2^^^  just  cancel  these  contributions  of  the  bound 

(2) 

states  the  energies  of  idiich  are  lower  than  E^'  . 6n  ^12)  becomes 


6n 


■ /' 


^ (:^,k2)  (3^,k2)  + ^ (x2,k2‘^)  (2.33) 

j 


This  leads  to  the  asymptotic  behavior: 


-9h  H X.. 

l6n  (Xj)l  < (const.)  e ^"2  2 as  x„->- 


(2.34) 
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where  h2  is  the  distance  from  the  real  axis  to  the  wave  vectors  correspon- 
ding to  the  highest  occi^ied  state  in  the  plane. 

Case  3:  above  the  energy  of  the  first  unoca;5)ied  state: 

the  bound  state  contributions  of  Eq.  (2.30)  then  cancel  the  residues  of  the 
(2) 

zeroes  of  1^2  ' vdiich  are  corresponding  to  the  occupied  bound  states.  6n 
becomes 


j 

c 1^22 

I 


'^2  '^2(x2.i^) 


1^2  (X2,k2^  ) ^2  (x2>^2^  ) 


where  j ' runs  over  all  poles  of  J^2^  ^ corresponding  to  the  unoccupied  bound 

(2) 

states  the  energies  of  vMch  are  lower  than  . This  leads  to  the  asynptotic 

behavior: 


|6n(3^)|  < (const.)  e''‘'^2  ^2  as  X2 


vdiere  h is  the  distance  from  the  real  axis  to  the  wavevectors  corresponding 
2 

to  the  first  unoccupied  st£ 


C.  Slnple  Interface: 

As  the  laigth  of  the  potential  (x)  goes  to  infinity  (N  -wo  ) , the 
system  shown  in  Fig.  2 looks  like  tvro  one-dimensional  periodic  lattices  joined 

at  the  origin  (Fig.  6).  In  this  case,  ^ analytic  function  of  k2 

(T) 

except  at  branch  points  vdiich  are  the  same  as  those  of  at  wavevectors 

corresponding  to  the  band  edges  of  a perfect  lattice  v/ith  the  potential  (x) 
and  at  the  infinitely  many  pole  singularities  vMch  occur  at  the  zeroes  of  ^22^^^ 
(k2)  (which  corresponds  to  the  existence  of  the  bound  states) . 

Because  of  this  analyticity  and  the  periodic  properties  of  the  integrand  of 
Eq.  (2.30),  we  can  shift  the  path  of  integration  from  the  real  axis  to  the  contour 
C as  indicated  in  Fig.  8.  The  change  of  the  electron  charge  density  at  X2  is 
(Appendix  B) 


(V . j 


^2  (x2,k2)  <1^2  (^{2,^2)  ) 


^^2'^'  j-  j 


(2.37 


The  second  term  on  the  right  hand  side  of  Eq.  (2.37)  is  the  contribution  from  all 

the  occupied  bound  states  which  are  located  in  the  gap  between  the  first  vmfilled 

and  last  filled  band  of  a perfect  lattice  with  potential  V2(x).  The  third  term 

(2) 

is  the  contribution  from  the  poles  from  the  zeroes  of  M22  ' \diich  are  located 
between  the  real  k2  axis  and  the  contour  C . 


Again,  we  can  discuss  Eq.  (2.37)  in  three  cases  as  we  did  for  the  isolated 


defect. 


Case  1:  lies  between  the  energies  of  the  first  unoccupied  state 

and  the  last  occupied  state:  the  second  and  the  third  terns  on  the  ri^t  hand 
side  of  Eq.  (2.37)  then  cancel  each  other  identically.  6n  (x2)  becomes 

f 

/__  \ _ I Jl.  ..  / _ U \.I.  U \ i 


f ^^2 

(X2)  = /dk2  *^2 

{ ^1^'’ 


(2.38) 


We  then  have  the  asynptotic  behavior: 


|6n  (X2)  1 < (const.)  e”^^  ^ as 


X2  ® 


(2.39) 


Case  2:  lies  below  the  energy  of  the  last  occupied  state:  tlic 

third  term  on  the  right  hand  side  of  Eq.  (2.37)  then  just  cancels  the  contribu- 

(2) 

tion  of  the  bound  states  the  energies  of  which  are  lower  than  E^  . iSn  6C2) 


becomes 


CX2)  ^ '^2  <X2*^)’i'2  ^2’*^^  (X2,k2-^)  (2.40) 

^ ^2^'^  j 

E^^^  (k2j)>EQ^^\ 


Wle  then  have  the  asymptotic  behavior: 


-2h^  30. 

|6n  (3^)1  < (const.)  e 2 ^ , as  35^ 


(2.41) 


Case  3:  lies  ^uXJve  the  energy  of  the  first  unoccupied  state:  the 

second  term  on  the  right  hand  side  of  Eq.  (2.37)  then  cancels  all  the  residues 
from  the  poles  corresponding  to  the  occupied  bound  states,  6n  (X2)  becomes 

f «12®  . r«L2® 

6n  (X2)  - /dk2  — 1^2  U2f^H2  (x  *^2)  -2-^1  ^res  — ij 

i ”22^^  j' 


^2  (X2.k2'^  ) 


/X  ’ 

\ * • 


(2.42) 


vhere  j runs  over  all  the  i'ol 


sponding  to  tlie  unocevpied  bound  states 


the  energies  of  vjhich  are  lower  ti’an  E^'  , We  then  have  the  asynptotic 

behavior: 


(3^)1  < (const.)  e‘^^2^  ^ as  3^ 


(2.43) 


We,  therefore,  have  concluded  that  the  electron  density  distribution 
moving  away  from  the  interface  approaches  exponentially  the  bulk  values.  The 
decay  constant  depends  ijpon  the  relative  position  of  the  band  structure  of  the 
perfect  lattices. 

III.  Mjre  General  Systems 

In  this  Section,  we  extend  the  theory  developed  in  Section  II  to  study  the 
following  question:  if  a point  perturbing  potential  is  introduced  into  a system 
\diLch  contains  a succession  of  n arbitrary  different  insulators  as  shown 
in  Fig,  9,  how  does  the  electron  density  perturbation  produced  by  this  distur- 
bance appear  across  successive  interfaces? 


./*■  >,.  j;--.  ■■-.S:.--V^TS»„  y-y,.  . ^Yp.:  ;V-  ^ J 


fM 


First,  let  us  consider  the  sisplest  case:  a system  with  just  one  interface. 

Vfe  take  the  perturbing  potential  (defect)  at  point  P a distance  to  the 
left  of  the  interface  as  indicated  in  Fig.  10. 

To  obtain  the  transfer  matrix  of  the  system,  we  will  follow  the  Saxon's 
method  as  we  did  before.  In  this  method,  we  only  replace  the  scattering  matrix 
at  X = -x^  (see  Eq.  (2.23)  ) by  the  defect  scattering  matrix 
Ihen  the  transfer  matrix  of  the  system  with  the  defect  is 

. r(2>  (s^(2)  r(2))^  (s^(2)  ^2)^  ^^(T>  RCa)/”^'!  „ „ 


and  that  without  the  defect  is 


M<«  . r(2)  (s^(2)  (s^(2)r(2))  (s^(2)  r(2)/-‘\-1  3_^<»  . (3.2) 


■Hie  analytic  properties  of  the  Mj.  yi^'^  similar  to  those  of  the 

(2) 

1^2  • both  are  analytic  functions  except  at  branch  points  and  poles  as  dis- 

cussed above  in  the  case  of  the  simple  interface.  However,  though  the  branch 

points  occur  at  the  same  k2  values  for  both  functions,  the  poles  occur  at 

(2) 

different  ^2  valvies  (the  zeroes  of  Mj-  22  ^ general  differ  from  those 

of  ) . 

Using  the  same  procedure  as  we  followed  for  the  case  of  the  sinrple  interface, 
we  have  the  following  form  ror  the  change  of  the  electron  density  distribution 
at  a distance  X2  to  the  right  of  the  interface  due  to  a defect  at  a distance 
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I 


to  the  left  of  the  interface 


6n  (xj^,  3^)  = - n (j^,  3^) 


J *2 


“r  ”12'^^ 

( — ' ~ - — ) ^\)2  (.^.^2)^2 


+ (3^,k2^)-^n®  (3^,k2j) 


2TTi  < ^ res 

J’ 


^,22^ 


’J^2  )i|>2(x2.k2'^  ) 


A • 

^2 


.kj^')  I 


(3.3) 


viiiere  the  subscript  I denotes  the  system  with  a defect;  the  notation  otlierwise 
being  the  same  as  in  Section  II. 

By  the  same  reasoning  vdxLch  we  used  to  discuss  Eq.  (2.37),  we  find  that 
n (Xj^,X2)  has  the  following  asynptotic  behavior-. 

|6n  (3£j^,3C2)|  < (const. )^  e 22^2^2  as  3^  ” (3.4) 

where  the  subscript  indicates  k^t 'constant  in  this  equation.  02 
is  the  decay  constant  in  the  insulator  2;  it  has  different  values  for  the 
following  three  cases: 


(T) 

Case  1:  When  the  energy  of  the  first  unoccupied  state  is  higher  than  ' 
and  the  energy  of  the  last  oca:5)ied  state  is  lower  than  , 

0.2  « h2°  . (3.5) 

(2) 

Case  2:  When  the  energy  of  the  last  occvqsied  state  is  higher  than  , 

02  = h2“  . (3.6) 

(2) 

Case  3:  When  the  energy  of  the  first  unoccupied  state  is  lower  than  , 

02  * ^2^  (3.7) 

If  in  case  (2)  the  bound  state  from  the  defect  is  the  last  occupied  state 
or  if  in  case  (3)  it  is  the  first  unoccupied  state,  then  02  depends  upon  the 
strength  of  the  defect  potential.  Otherwise  02  depends  upon  the  potentials  of 
the  two  lattices  only. 

We  now  want  to  discuss  the  modification  of  6n  (xj^,X2)  as  the  defect  is 
moved  away  from  the  interface  vdiile  X2  is  fixed.  To  obtain  the  functional 
dependence  of  6n  (xj^,:^)  i^x>n  Xj^  , we  have  to  change  the  fundamental  inter- 
val in  Fig.  2 as  follows:  it  is  taken  to  be  a perfect  lattice  with  potential 
Vj^  (jO  with  extended  impurities  from  x = 0 to  x = N of  the  potential  V2  (x) 
as  indicated  in  Fig.  11.  The  potential  in  this  fundamental  interval  is  given  by 
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V(x) 


V|^(x) 

-d/2«<0 

V2  (x) 

0<  x<N 

Vj^  (X) 

N<  x<d/2 

(3.8) 

Using  the  same  prcx:edure  as  befcare,  we  obtain  the  transfer  matrix  of  this 
system  by  Saxon's  method; 


S (l).-l^(l)g  (1) 


(3.9) 


and 


(3.10) 


All  the  notation  is  the  same  as  before,  except  for  the  si^erscript  (1)  which 
denotes  k^  as  an  independent  variable. 

From  the  Appendix  D,  we  obtain  the  following  asymptotic  behavior  of  6n  (xj^.i^) 
with  fixed  X2 


“2ot  y 

[fin  (x^,X2)|  < (const.)  e 1"^  as  ^ 

*2 


(3.11) 


where  otj^  is  the  decay  constant  in  the  insulator  1.  It  is  given  by 


3 
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if  the  energy  of  the  first  unoccupied  state  is  hig’ner  than  E 

/I  \ o 

energy  of  the  last  occupied  state  is  lower  than  E ' ' ' 


“l  “ ^1 


if  tlie  energy  of  the  last  occupied  state  is  higher  than  E 


a,  = h 


if  the  energy  of  the  first  unoccupied  state  is  lower  than  E 


where  , the  branch  points  of  the  energy  surface  is  located 

between  the  first  unfilled  and  last  filled  band  of  a perfect  lattice  with  poten- 
tial (x).  h^°  , h^^  and  h^^  are  the  distances  from  the  real  axis  to 

the  wavevectors  corresponding  to  the  branch  points,  the  last  occupied  state, 
and  the  first  unoccupied  state,  respectively. 


From  Eqs:  (3.4)  and  (3.11),  we  have  the  following  asynptotic  behavior 
of  6n  (^.^): 

-2a,  -2^2^ 

1 6n  ( Xj^ , :^)  I < (const . ) e e as  Xj^  -*  <» 


(3.12) 


I 

I namely,  a system  with  n interfaces.  This  system  has  a succession  of  n + 1 

[ arbitrary  different  sublattices  as  showi  in  Fig.  9.  The  potential  in  the 


vth  sublattice  is  denoted  by  v^(x).  All  the  sublattices  except  for  those 
at  the  ends  are  of  finite  length,  denoted  by  x^;  the  end  lattices  are  infinite. 
We  put  an  isolated  defect  in  the  first  sublattice  at  a distance  x^  to  the 
left  of  the  interface  vhich  divides  the  first  and  second  sublattices.  We 


«ant  to  determine  the  change  of  the  electron  density  distribution,  6n  C^.X2  ••>^n+l)’ 
in  the  last  (v=nfl)  sublattice. 

If  the  length  of  each  sublattice  is  fixed,  we  have  from  the  above  discussion 
the  following  asymptotic  behavior  of  6n  .x^  . . . with  respect  to  x^^ 


|6n  (x^,X2,  ® ® 

as  x^ 

''nfl  " ” (3.13) 

We  now  want  to  find  the  functional  dependence  of  6n  0c^pc2, . . . upon 

x^  for  v“2,3...,n.  First,  we  fix  all  x^  except  a particular  one  y . In  this 
case  we  can  think  of  the  lattice  y as  a host  lattice  while  the  other  lattices 
are  thought  of  as  extended  iinpurities.  We  choose  the  fundamental  interval 
shown  in  Fig.  12  and  follow  the  same  procedures  as  we  followed  in  Sec.  II.  We 
have  (see  Appendix  D) 
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l<5n  I (const.)  x ^ x ^ 


X , X . e as  X -»oo 

Vri-1”“’  nf  1 ® y 


(3.14) 


where  a^,  the  decay  constant,  has  the  similar  definition  as  before. 

Then  we  have  the  following  results  for  the  change  of  the  electron  density 


distribution  due  to  a defect: 


n + 1 


|<Sn  (xj^,x2 x^l)l  < (const.)  j~| 


e as  X -»•  “ 


(3.15) 


v=  1 


where  a is  the  decay  ccmstant  in  the  sublattice  v . Its  value  is  given  by 


a , = h o , 

V V 


if  the  energy  of  the  last  occupied  state  is  lower  than  and  the  energy 

of  the  first  unoccupied  state  is  higher  than  E . 


V V 


if  the  energy  of  the  last  occupied  state  is  higher  than  ^ . 


a . = h 


if  the  energy  of  the  first  unoccupied  state  is  lower  than  E 


^Aiere  , the  branch  points  of  the  energy  surface  , is  located 

between  the  first  unfilled  and  last  filled  energy  band  of  a perfect  lattice 
with  potential  v (x) . h ° , h ^ , and  h ^ ar  e the  distances  froci 

V V V V 

the  real  axis  to  the  wavevectors  corresponding  to  the  branch  points,  the  last 
occupied  state  and  the  first  xiioccupied  state,  respectively. 

IV.  SunnHry 

The  perturbation  of  the  electron  density  distribution  due  to  an  isolated 
defect  appearing  across  nfl  successive  arbitrary  different  lattices  decays 
exponentially.  The  decay  constant  in  a particular  sublattice  p dopopHR 
on  tlie  position  of  the  branch  points  E^^^^  (vMch  are  located  between  the 
last  filled  and  first  unfilled  energy  bands  of  a perfect  lattice  with  potential 
v^  ( 3d . ) relative  to  the  energies  of  the  highest  occiqsied  state  and  the  last 
unocciq)ied  state  of  this  S3rstan. 

There  are  three  cases: 

Case  1:  If  the  energy  of  the  last  occupied  state  is  lower  than  E^^^^ 
and  the  energy  of  the  first  unoccupied  state  is  higher  than  E^^^\  the  decay 
constant  equals  the  distance  fran  the  real  ascis  to  the  wavevectors  corresponding 
to  the  branch  points  E^^^^  in  the  plane. 

Case  2:  If  the  energy  of  the  last  occiqjied  state  is  higher  than  E^^^^  , 


the  decay  constant  equals  the  distance  from  tlie  real  ascis  to  the  vavevectors 


corresponding  to  the  last  occupied  state  in  the  k^  plane. 

Case  3:  If  the  energy  of  the  first  unoccvipied  state  is  lower  than  E 


(y) 


the  decay  constant  equals  the  distance  from  the  real  axis  to  the  wavevectors 
corresponding  to  the  first  unocciqDied  state  in  the  k plane. 
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APPENDIX  A 


In  this  appendix,  we  want  to  discuss  the  analytic  properties  of  the  trans- 
fer matrix  M (see  Eq.  (2.18))  for  a system  with  an  interface.  (For  simplicity 
of  notation,  the  superscript  (2)  will  be  suppressed  in  this  appendix.) 

To  calculate  , we  have  to  find  a transformation  P such  that  G is 
in  diagonal  form: 


P = 


% (&+-o),  % (3-a) 

, R2ie‘^2 


(A.l) 


v^iere  a and  3 are  defined  in  Eq . (2 . 23)  . G can  be  written  as 


G = 


(A.  2) 


where  and  X2  are  defined  in  Eq.  (2.23)  and  vhere  P"^  is  the  inverse 
matrix  of  P given  by 

p"l  = -i-  I 

Here  |P|  is  the  determinant  of  P given  by 


-ik2 

(p|  - R21®  “• 


(A.  4) 


Prom  Eq.  (A.l)  and  Eq.  (A,2),M  can  be  written  as  Eq.  (2.22). 
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I 


We  now  \igant  tx)  discuss  the  anal3rtic  properties  of  ^ system 

with  an  Interface  (N>«>).  Fran  Eq.  (2.22),  given  by 


^2 

*4.1^12®^^  (Ai”  - A2^)  +\2^ 

(&+u)  A2^  - (B-a)Aj^^  ] 

^2 

^1^2®  ^^1^  ■ ^2^^  +^22^ 

(&+a)  A2”  - (e-a)  . 

Before  proceeding  with  the  discussion,  let  us  give  the  physical  meaning  of 
and  A2.  By  comparing  the  Saxon's  method  with  the  Heine's  method  we  have 
the  following  relations  between  and  k2: 

(A.  6) 

vhere  satisfies  the  following  relation: 

£^2)  (1^)  = (k^)  . (A.  7) 

Here  and  E^^^  (k2)  are  the  energy  surfaces  of  the  perfect  lattices 

with  the  potential  Vj^(x)  and  V2(x),  respectively. 


Following  those  relations  , A^^  and  A2  can  be  written  as: 


Let  us  now  return  to  the  discussion.  ^ analytic  function  of 

CD 

k2  , excqjt  at  its  branch  points  vMch  are  the  same  as  those  of  E'  ' (k2)  and 
\diich  occur  for  those  valvies  of  '^ch  correspond  to  being  mT;n  = 0, 

+1,  +2,  — (a*0),  and  has  infinitely  many  poles  \diich  occur  for  those  valxies 
of  k2  v^ch  correspond  to  the  real  , i.e.  those  poles  occur  at  the  points 
corresponding  to  the  points  in  the  allowed  bands  of  and  the  branch 

points  occur  at  the  points  corresponding  to  the  band  edges  of  E^^^  (k^^) . 


APPENDIX  B 


Me  consider  here  a one-dinensional  lattice  to  xdiich  extended  iiipjritles 
have  been  added  (see  Fig.  2) . The  eigenstates  of  the  system  are  the  scatterizig 
states  ^ (x)  and  the  bound  states  (Ji®  (x,k2^).  Here  (+)  and  (-) 

denote  states  with  outgoing  and  incoming  scattered  waves,  respectively.  Out- 
side the  iiqpurity  region,  the  scattering  states  are  appropriate  linear  combination 
of  Bloch  waves 


(x) 


^ I 4-2(x.k2)  x>o 


(x)  = 

^ (Moo 


(B.l) 


i|.2(x  ,-k2) 


x<-  N 


■Hie  coefficients  and  ^2^^^  elements  of  the  transfer  matrix. 

At  point  X2  (X2>0)  , the  change  of  electron  density  due  to  those  impuri- 


ties is 


6n(x2)=  nj(x2)  - n(x2) 


dk2  '^-k2^"^  (X2,k2^) 


T 


I 

b 


I 


.JL 

b 

H 

b 


/ 


dk2 


‘^2 


(2) 


TT 

'S’ 


^2 


(2) 


'^2^^’^^  ■'■  (X2.k2'^) 


(B.2) 


(2')  (2) 

Fran  ^jpendix  A,  we  know  11^2  '^^2  ^ analytic  function  of  k2,  except 

(2) 

at  an  infinite  number  of  poles  vdiich  occur  at  the  zeroes  of  three 

sets  of  the  branch  points  which  occur  at  the  points  k2°  , k2'^  and  k2^  , corres- 

ponding  to  the  branch  points  of  the  energy  surface  (k2)  , and  at  the  upper 
and  lower  band  edges  of  a perfect  lattice  with  the  potential  Vj^(x) . If  we 
deform  the  path  of  integration  of  the  integral  in  Eq.  (B.2)  from  the  real  axis 
to  the  contour  C just  below  the  branch  points  k2°  , the  integral  in  Eq.  (B.2) 
can  be  written  as  a single  integral  along  the  contour  C plus  the  residues  at 
the  zeroes  of  M22'  ' which  are  located  between  the  real  axis  and  the  contour 
C . This  results  because  of  the  calculation  of  the  vertical  portions  of  the 
contour  (this  is  due  to  the  periodicity  of  the  integrand  in  k2) . 6n(x2)  becomes 

(2) 


■/ 


6n(x2)  = I dk2 


^2' 


M22 


(2) 


(X2.k2'^) 


-2Tri  res 

j’ 


^2 


(2) 


^2 


jjy  \p2  (5^2’^^  ^ "^2  ^ 


(B.3) 
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Again,  from  the  appendix  A , vge  know  the  integral  in  Eq.  (B.3)  has  the 
following  asynptotic  behavior: 

I /dk2  ^2)  "^2  I ^ (const.)  e ^ as 

•/  Moo 


(B.4) 


vhere  h2°  is  the  distance  from  the  real  axis  to  the  branch  points  k2°. 

The  bound  state  contributions  to  6n(x2)  are  denoted  by  Ej  n®(x2,k2^), 

vhere  j runs  over  all  the  oca.?)ied  bound  states  and  n (X2,k2'^)  is  bounded 

• • 

by  (const.)  e:q)(-2h2'^X2) . Here  h2^  is  the  distance  from  the  reax  axis  to  the 

wave  vectors  k2-^  corresponding  to  the  bound  state  j with  energy  = E^^\k2^). 

The  residues  at  the  zeroes  of  ^2^^^  denoted  by  the  second  svnmation 

on  the  right  hand  side  of  Eq.  (B.3).  j runs  over  all  the  zeroes  of  ^2^^^ 

vhich  are  located  between  the  real  axis  and  the  contour  C . The  residue  at 

the  zero  of  ^ is  bound  by  (const.)  e3q)(-2h2'^  X2),  where  h2''  is 

the  distance  from  the  real  axis  to  the  k2'^  . 

We  will  discuss  Eq.  (B.3)  in  the  following  three  situations: 

(1)  if  the  energy  of  the  last  occupied  state  is  lower  than  and 

the  energy  of  the  first  unoccupied  state  is  higher  than  E^^^^  , 

j and  j'  on  the  right  hand  side  of  Eq.  (B.3)  run  over  the  same  set  of  k2  . 

The  bound  state  contributes  to  precisely  cancel  the  residues  at  the  zeroes 

(2)  1 

of  1^2  • Therefore,  6n(x2)  has  the  following  ^lsyIIptotic  behavior: 


-2h2°X2 

|6n(x2)  |<(const.)e  as  X2-«» 


(B.5) 


If  the  energy  of  the  last  ocojpied  state  is  higher  than  E , the 
(2) 

residens  at  the  zeroes  of  1^2^  ^ between  the  real  axis  and  the  contour  C will 
cancel  the  contributions  of  the  bound  states  with  energy  levels  below  . 

Eq.  (B.3)  becomes 

6n(x2)  = y dk2 )J^2  ^'^>^'^2  ^ ” (X2.k2^) 

C ^2  j 

(B.6 


In  the  forbidden  region  vAiere  the  energy  is  higher  than  E^^  ' , we  know 
from  the  Kramer's  plot  that  the  higher  the  energj/of  a bound  state,  the  smaller 


the  distance  from  the  real  axis  to  the  wave\iectors  corresponding  to  this  bound 
state.  Then  6n  (X2)  has  the  following  asynptotic  behavior: 


1 6n  (X2)  |<  (const.)  e 


-2h2X2 


as  Xo->«> 


vdiere  h2  is  the  distance  from  the  real  axis  to  the  vravevectors  corresponding 
to  the  highest  occupied  state. 

(2) 

(3)  If  the  energy  of  the  lowest  unocctqjied  state  is  lower  than  E^'  , 

the  second  term  on  the  right  hand  side  of  Eq.  (B.3)  (the  contributions  from  the 
occupied  bound  state)  will  cancel  the  residues  at  the  poles  corresponding  to 
the  occupied  bound  states.  6n(x2)  becomes 


,r(x2)  - / dk2  '<^2 


r V 

“2'n'i  ^ res  *1^  (X2fk2'^  ) 

j.  L“22  J 


(B.8) 


\Aiere  j'  runs  over  all  the  zeroes  of  1^2^  ^ corresponding  to  the  unoccupied 

(2) 

bound  states  with  energy  levels  below  E ' . 


In  the  forbidden  region  where  the  energy  is  lower  than  £ 


we  know 


from  the  Kramer's  plot  that  the  higher  the  energy  of  a bound  state,  the  larger 
the  distance  from  the  real  axis  to  the  wavevectors  corresponding  to  this  bound 
state.  Then  6n(x2)  has  the  following 


|6n(x2)l 


-2h2Hc2 

< (const.)  e ^ 


as  Xo 


(B.9) 


v^iere  h2  is  the  distance  from  the  real  axis  to  the  wavevectors  corresponding 
to  the  lowest  unocci;q)ied  state. 


I 


APPEM)IX  C 


In  this  appendix,  we  want  to  derive  Eq.  (3.11) . If  we  define  A,  A^.  and 
6 as  follows: 


the  lij  and  M defined  in  Eqs.  (3.9)  and  (3.10)  can  be  written  as: 
ri,  = B . A, 

M = B • A (C.2) 

Again,  for  sinplicity  of  notation,  we  suppress  the  superscript  (1)  in  this  appen- 
dix. 

The  change  of  the  electron  density  distribution  at  the  point  ^2  due  to  the 
defect  at  the  point  is 

IT 

ir  1 1 i -< 

'F 

+ (X2,ki-^)  -^n®(x2,kj^-^)  , (C.3) 

j j 
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1 


;*ere  (x)  are  the  scattering  states  of  the  system  with  the  defect; 

(+) 


(X2)  =■ 

1 

<»I,22)^ 

f^22  *^1  ^^^’^^^”^21  J 

(-) 

(X2) 

1 

(Mj-  22^^ 

K.12  -^,22 

(+) 

'^■+kj^  (x)  are 

the  scattering  states  of  the  system  without  the  defect: 

(X2)  = 

(M22)^ 

(C.5) 


K2  V + A22'i'j^(x2.-lcj^)  ] 

Hie  integrand  of  the  Eq.  (C.3)  I (x^,X2,k^)  can  be  written  as: 

Kxj.xj.k^)  = [b22»i  (x2.ki)-B2i»i(x2,-k^)  j 

-Aj  22'I'i(>=2’-V  ] - -^  [ 

" (^2^  22  **1,12*22A2*i  22'  [°22*‘l*’^''U'‘®21*l^’^’‘’^'  ] 


*’'12® 


2ik^Xj^ 


^22 


N-x^  -ik^(x2+l) 


N-: 


ik^Cx^+l) 

e (x2.-k^) 
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X2+1  2ilCj^Xj^ 


+T22G12 


X2+I) 


(T,,G, 


X2+I  2ikj^x^ 


12''21 


-HT22G22""^'>1  1^21  ^12 


N-X2  X2+I  N-X2  X2+I 


+G22  G22 


vAiere  GeS^^R. 

I(x^,X2.1tj^)  lias  the  same  analytic  praperties  as  the  integrand  of  Eq.  (B.2). 
We  can  deform  the  path  of  integration  of  the  integral  in  Eq.  (C.3)  to  a contour 
C just  below  the  branch  points  of  E^^^  (kj^)  as  we  did  in  Appendix  B.  We  obtain 
the  asynptotic  behavior  of  the  following  integral  along  the  contour  C with  fixed 


I j dkj^I(Xj^,X2,k^)| 


-2h,°x, 

< (const.)  e as  X, 


(C.7) 


where  h^°  is  the  distance  from  the  real  axis  to  the  wavevectors  corresponding  to 
the  branch  points  of  E^^^  (kj^)  . 

Using  Eq.  (C.7)  and  following  the  same  reasoning  as  we  did  in  Appendix  B, 
we  can  obtain  Eq.  (3.11). 


i 


APPENDIX  D 


In  tMs  appaidix,  we  want  to  derive  Eq.  (3.14).  For  sinplicity,  we  consider 
here  the  system  with  two  interfaces  as  shown  in  Fig.  9 instead  of  n inter- 
faces. There  are  three  sublattices  1,  2 and  3 with  the  potential  v^(x), 

V2(x)  and  v^Cx) , respectively.  Since  vje  want  to  obtain  the  asymptotic  behavior 
of  the  change  of  the  electron  density  distribution  in  the  sublattice  3 due  to 
a defect  in  the  sublattice  1 as  the  length  of  the  sublattice  2 X2  approaches 
to  infinity,  we  can  tliink  of  the  sublattice  2 as  a host  lattice  and  k2  as 
an  independent  variable  while  the  other  sublattices  1 and  3 are  thought  of 
as  extended  inpjrities.  The  fundamental  interval  of  this  system  is  shown  in 
Fig.  11.  The  lengtlies  of  the  sublattices  1 and  3 are  denoted  by  and 
respectively. 

The  transfer  matrices  for  this  system  with  and  without  the  defect  are  given 


M = B • A 


(D.l) 


vdiere 


B = S, 


(SiK) 
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7 


-1  X-3  -1  t 


A . s ' (S,K)  3 S.^  (S^R)  3 


(D.2) 


Here  R arid  K are  the  scattering  matrices  due  to  a localized  inpurity  at  the 


origin  in  a perfect  lattice  2 with  the  inpurity  strength  v^(x)  and  v^Cx) , 
respectively  and  \diere  the  point  x^  is  in  the  sublattice  3 and  the  point  x^^ 


is  in  the  sublattice  1. 


The  change  of  the  electron  density  distribution  at  the  point  x^  due  to 
the  defect  at  the  point  x^  is  given  by 

IT 

F 

6n(x^.X2.X3)  = j dk^  ^ ^\,12^2~h2^.22^ 

•'tt  x,22z2 

F 

+ ^ n^®  (X2,k2'^)  - (X2,k2-^)  < 


\diere 


lR^(2^,+2x.)  r 

^,l2^2~hA,22  = ® Kl^r\l^21^ 


^12^  ^12 


^'^11^2'^'^12^r\l^22'\2^21^  ^12  ^22 


(Ti2R22~^2'^22^  ^22  ^’22  ^ , 


1 


ij-Xj  ik2>b 


and  \diere 


£^-x-  r ik„  x-  x, -1  2ilc,x«  x^  x--l 

••■l,22»22  = °21  ”21  I®  <”ll  <=11  =^^■^”12  ”21  > 

£,-x,  "ik«  X-  x,-l  2ik«x,  x,  x,-l 

(T11G12  +T^2”22  <”11  ”11  ' -*12  ”22  XT2A2 

T g V’h.  )|  (D,  \ 21^2=^  >=3  x^-1 

■^22^22  “'J  L®  *^11  6 ^12  ^21  ^ ^\l^l2 


£l-Xi 


■^2^22 


£,-x^  -ik,  X,  ^.1  2ik2^  X3  x^-1 

/ e ^ (Dl,  ^1  e -m,.  G..  ) 


'12  ^22  ^ ^^A2 


^,-x,  I ik^(2£,+2£p) 

^2”22  ^ ^ 


(D.4) 


As  and  go  to  infinity,  the  integrand  of  the  integral  in  Eq.  (D.3), 

I (x^,X2,X2,k2)  has  analytic  properties  similar  to  those  of  the  integrand  of 
the  integral  in  Eq.  (B,2).  It  is  an  analytic  function  of  k2  , except  at  the 
branch  points  corresponding  to  the  branch  points  of  E^^^(k2)  and  the  band  edges 
of  the  perfect  lattices  vdth  the  potential  Vj^(x)  and  V2(x)  and  at  the  poles 
v^ch  occur  at  the  zeroes  of  1^2  ^ 22‘  deform  the  path  of  inte- 

gration from  the  real  axis  to  a contour  C just  below  the  branch  points  of 

(2) 

E (k2) . We  obtain  this  asynptotic  the  contour  C for  fixed  x^^  and  x^ : 


-40- 


X 


iydk2l(Xj^.X2,X3.k2)|<  (const. 


-2h2°X2 
e as 


(D.5) 


\(*iere  h2°  is  the  distance  from  the  real  axis  to  the  branch  points  of  (k2) 
Again,  using  Eq.  (D.5)  and  following  the  same  reasoning  as  we  did  in 
^jpendix  B,  we  can  obtain  Eq.  (3.14) . 
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Integration  contour,  C,  for  the  integral  in  Eq.  (2.40).  k2 
is  a branch  point  of  e(2)  (k2)  in  the  k2  plane;  k2“  ana 
wavevectors  corresponding  to  the  band  edges  of  a perfect  lattice 
with  the  potential  vj^Cx).  The  tvx)  closely  spaced  seqiaences  of 
points  indicate  the  positions  of  wavevectors  associated  with  the 
bound  states  due  to  the  N cells  with  potential  v-^,  embedded  in 
the  infinite  potential  Vry. 
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